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Recent data from p+p and p+Pb collisions at the Large Hadron Collider (LHC), and d+Au 
and ^He+Au collisions at the Relativistic Heavy Ion Collider (RHIC) reveal patterns that—when 
observed in the collision of heavy nuclei—are commonly interpreted as indicators of a locally equi¬ 
librated system in collective motion. The comparison of these data sets, including the forthcoming 
results from p-|-Au and p-l-Al collisions at RHIC, will help to elucidate the geometric dependence of 
such patterns. It has recently been shown that A-Multi-Phase-Transport-Model (AMPT) can de¬ 
scribe some of these features in LHC data with a parton-parton scattering cross section comparable 
to that required to describe A + A data. In this paper, we extend these studies by incorporating 
a full wave function description of the ®He nucleus to calculate elliptical and triangular anisotropy 
moments V 2 and vs for p-|-Au, d-|-Au and ^He-I-Au collisions at the RHIC top energy of 200 GeV. 

We hnd reasonable agreement with the measured V 2 in d-|-Au and ^He-|-Au and vs in ®He-|-Au for 
transverse momentum (pr) fS 1 GeV/c, but underestimate these measurements for higher values of 
Pt- We predict a pattern of coefficients {v2, wa) for p-|-Au, dominated by differences in the number 
of induced local hot spots (i.e. one, two, or three) arising from intrinsic geometry. Additionally, 
we examine how this substantial azimuthal anisotropy accrues during each individual evolutionary 
phase of the collision in the AMPT model. The possibility of a simultaneous description of RHIG- 
and LHG-energy data, the suite of different geometries, and high multiplicity p+p data is an exciting 
possibility for understanding the underlying physics in these systems. 

PACS numbers: 25.75.Gz, 25.75.Gz.Ld 


I. INTRODUCTION 

Recent analyses of data from p+p and p-|-Pb colli¬ 
sions at the LHC, as well as from d-|-Au and ^He-I-Au 
at RHIC have revealed the existence of azimuthal parti¬ 
cle correlations reminiscent of those observed in A -|- A 
collisions [2J HI O [T31 HH US] . In the latter, this signal 
has been interpreted as evidence for the liquid-like na¬ 
ture of the Quark-Cluon Plasma (QCP)—a locally equi¬ 
librated, strongly coupled medium undergoing hydrody¬ 
namic expansion with viscosity near the quantum lower 
bound. However, these results were largely unexpected in 
p(d)+A, long considered control systems to understand 
initial-state effects in heavier systems. 

The success of nearly inviscid hydrodynamics in de¬ 
scribing A -|- A bulk observables makes it natural to ask 
if droplets of QCP are being formed in small systems, and 
if the created medium is sufficiently long lived to equi¬ 
librate locally and translate initial spatial anisotropies 
into final-state particle momentum correlations. How¬ 
ever, this is not the only possibility, with potentially dif¬ 
ferent physics being able to account for these measure¬ 
ments mm [T8| . Further insight into this matter will 
come from the confrontation of different model calcula¬ 
tions with the full data sets available. 

It has been found that azimuthal two-particle corre¬ 
lations from high multiplicity p+p, p-l-Pb, d-|-Au, and 
^He-|-Au events exhibit an enhancement around « 0 
(i.e. near-side), even when the particles have a large sep¬ 
aration in pseudo-rapidity {Arj > 2), where jet contribu¬ 
tions are expected to be minimal. There is an additional 


enhancement from p+p and peripheral p-|-Pb (d-|-Au) to 
central p-|-Pb (c?-|-Au) around A(j) « tt (i.e. away-side) 
that has been interpreted as the full azimuthal continu¬ 
ation of elliptical and triangular flow coefficients, V 2 and 
V 3 . Alternative readings of the away-side pattern include 
modification of the dijet correlations for the most central 
d-l-Au collisions [3]. 

Predictions for LHC-energy p-|-Pb anisotropies using 
nearly inviscid hydrodynamics provide a reasonable 
description of the flow coefficients measured at the LHC. 
However, as expected, an exact quantitative description 
depends on on the shear viscosity to entropy density ra¬ 
tio (rj/s) and the details of how initial geometry is mod¬ 
eled, for which there are quite different possibilities in 
p + A collisions—see for example Ref. [HI US] . There 
are also competing calculations where final-state QCP or 
flow effects are deemed negligible, and it is initial-state 
glasma diagrams that give rise to the correlations Ba¬ 
in d-l-Au and ^He-|-Au collisions, the initial geometry is 
dominated by the spatial separation of the two nucleons 
in the deuteron, reducing differences between models of 
geometry. For this case, nearly inviscid hydrodynamic 
calculations give a reasonable description of the experi¬ 
mentally extracted flow coefficients [ZlllllD]- 

However, questions regarding the validity of the near- 
inviscid hydrodynamic calculations have been raised in 
terms of the expansion around steep energy density 
gradients in these small systems mi unis]- It is 
thus quite interesting that incoherent elastic parton- 
parton scattering—as implemented in A-Multi-Phase- 
Transport-Model (AMPT) [T3] —can effectively repro- 
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duce the long-range azimuthal correlations m and V 2 co¬ 
efficients [ 8 ] observed in high multiplicity p+p and p-|-Pb 
events at the LHC. Notice, however, that in the case of 
V 2 in p-|-Pb, a good reproduction of the measured values 
is only achieved for px ^ 2 GeV/c, above which the cal¬ 
culations underestimate the data. These AMPT results 
were obtained using a parton scattering cross section of 
cr = 1.5 — 3.0 mb, and incorporating the so-called string 
melting mechanism in the model (thus including a time 
stage dominated by parton-parton scattering). 

These results raise the question of whether a simi¬ 
lar description can be achieved for different collision ge¬ 
ometries at the RHIC energy scale. In particular, we 
use AMPT to simulate p-l-Au, d-l-Au and ^He-I-Au at 
y/s^ = 200 GeV since they have been proposed as an 
excellent testing ground to disentangle the properties of 
the medium created in small collision systems [20) . In 
this paper, we begin by describing the AMPT model and 
the methodology used to compute azimuthal anisotropies 
of final state hadrons with respect to the participant 
plane. We then compare our V 2 and results in d-|-Au 
and ^He-I-Au to available data, and present predictions 
for V 2 and in p-|-Au. Finally, we discuss our results 
and provide some conclusions. 


II. METHODS 

The AMPT event generator m is a useful tool for 
the study of heavy-ion collision dynamics. The AMPT 
model uses the HIJING model m just to generate 
initial conditions via Monte Garlo Glauber, Zhang’s 
Parton Gascade (ZPG) to model partonic scattering, 
and A-Relativistic-Transport (ART) to model late stage 
hadronic scattering. We utilize the AMPT model with 
string melting turned on, such that a stage with parton- 
parton scattering is included and subsequent hadroniza- 
tion is described via a coalescence model. In this coa¬ 
lescence model, quark-antiquark pairs and sets of three 
(anti) quarks in close spatial proximity are grouped to 
form mesons and baryons, respectively. 

In order to better understand how AMPT translates 
anisotropies in the initial geometry to anisotropies in 
final-state momentum, we modified the internal Monte 
Garlo Glauber to more closely resemble the standard ap¬ 
proach described in [T71 HH] , and used in Ref. |20] ■ The 
position of nucleons in each colliding nucleus is sampled 
from the appropriate wave function on an event-by-event 
basis, after which a nucleon-nucleon inelastic cross sec¬ 
tion of 42 mb is used to geometrically determine which 
nucleons were wounded in the collision. In the case of the 
deuteron, coordinates are sampled from the two-nucleon 
Hulthen wavefunction; in the case of ^He, coordinates 
are obtained with Green’s function Monte Garlo calcula¬ 
tions using the AVI 8 -I-UIX model of three-body interac¬ 
tions [TO] . 

We have run approximately 10 million central (i.e. im¬ 
pact parameter b < 2 fm) AMPT events with a parton- 


parton scattering cross section u = 1.5 mb for each sys¬ 
tem at y/s~^ = 200 GeV. There are numerous publi¬ 
cations utilizing the AMPT model to describe A + A 
collisions at RHIG, quoting a full range of input cross 
section values ranging from a = 1.5 mb up to cr = 10 
mb mm [29]. For this study, we have chosen the 
smallest value from this range to understand how a min¬ 
imum parton scattering stage contributes to collective 
motion in these small systems. Table |l| summarizes the 
mean number of nucleon participants for each system in 
AMPT, as well as the corresponding yield of partons at 
the end of the parton scattering stage and the yield of 
final-state hadrons at the end of the hadronic scatter¬ 
ing stage. Note that all partons reported by AMPT at 
the end of the parton scattering stage are quarks and 
anti-quarks, and no gluons, which are then input to the 
particular coalescence calculation for hadronization. 

TABLE I. Particle production and eccentricity in central 
AMPT small-system collisions. For each collision system, 
we show the mean number of participant nucleons per event, 
the mean nnmber of partons (i.e. quarks and antiquarks) 
at freeze out, the mean number of hadrons after the hadron 
cascade, and the mean elliptical and triangular initial state 
eccentricities. 


System 

(y^part) 

partons) 

(y^hadrons ) 

(£2) 

(£3) 

p-l-Au 

10.45 

182 

131 

0.24 

0.16 

d-l-Au 

18.3 

336 

233 

0.57 

0.17 

®He-bAu 

22.3 

446 

326 

0.48 

0.23 


For each collision system, we compute the eccentricity 
and participant plane angle '^n on an event-by-event 
basis using the initial-state coordinates (r^, (pi) of partici¬ 
pant nucleons with the spatial distribution of a Gaussian 
of width a = 0.4 fm, as follows for n = 2, 3 

yj (r^ cos(n(/)))^ -I- (r^ sin(n(/)))2 

^n. = - ~ -- , ( 1 ) 

(^ 2 ) 

^ ^ atan 2 ((r^ sin(rt 0 )), (r^ cosinp))) ^ ^ 
n n 

Average values of £2 and £3 for central p-|-Au, d-|-Au, and 
^He-|-Au are shown in Table [H 

Having measured from the initial-state geome¬ 
try, we compute the second and third order azimuthal 
anisotropy moments V 2 and U 3 of final-state unidentified 
charged hadrons within \ri\ < 2 , with respect to the par¬ 
ticipant planes, as follows 

= (cos[n((p - 4'„)]). (3) 

In addition to extracting anisotropy moments with re¬ 
spect to the participant plane, we are also interested in 
qualitatively examining the long-range azimuthal corre¬ 
lations of these hadrons. To that end, we follow an anal¬ 
ysis procedure similar to that put forth by the ATLAS 
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experiment for p+Pb collisions in Ref. [T3]. We take all 
final-state charged particles and consider all pairs sepa¬ 
rated by 2.0 < |A? 7 | < 3.0 within a common px bin. We 
then form a long-range two-particle azimuthal correlation 
function for each px bin: 


C(A<^,Pt) 


1 dN{pT) 
Ntrig dA(j) 


( 4 ) 


Shown in the top panel of Figure |T] are the two-particle 
correlations from p+p and central ^e-|-Au collisions for 
particles within 0.9 < pr < 1.04 GeV/c. As expected, a 
flat near-side (around |A0| « 0) distribution is observed 
in p+p since two particles from the same jet fragmen¬ 
tation or resonance decay are very unlikely to be sepa¬ 
rated by more than two units of pseudo-rapidity. There 
is also a significant p+p away-side enhancement (around 
|A^| « tt) from jets back-to-back at leading order in 
azimuth and with a significant pseudo-rapidity swing be¬ 
tween them from an imbalance in the initial parton mo¬ 
mentum fractions Xi and X 2 - 

In contrast, in the ^He-|-Au distribution we observe a 
prominent near-side peak and a nearly two-fold enhance¬ 
ment of the away-side yield relative to p+p. Assuming 
the jet contribution to be the same in p+p and ^He-|-Au, 
we subtract the two distributions, thus obtaining the re¬ 
sult shown in the bottom panel of Figure We choose 
to use p+p data for the jet subtraction, analogously to 
what has been done in experimental data with peripheral 
p-l-Pb (d-l-Au) events [U [3 [13] . It is important to note 
that we do not use long-range two-particle correlations 
to extract anisotropy moments, but rather to provide a 
qualitative comparison with experimental results where 
similar correlations are presented. Hence, all results pre¬ 
sented in the remainder of this article are computed with 
respect to the participant plane using initial state geom¬ 
etry, as previously described. 


III. RESULTS 

The resulting anisotropy moments V 2 and V 3 for 
AMPT d-|-Au central collisions, computed with respect 
to the participant plane as described in Section II, as a 
function of px are shown in Figure We observe a sub¬ 
stantial V 2 that rises with px to around 10% at px ~ 1-0 
GeV/c, after which it levels off and even exhibits a slight 
decrease. The V 3 coefficients exhibit a similar px depen¬ 
dent trend, though substantially smaller values. Shown 
for comparison are elliptic flow measurements using the 
participant plane method for central (0-5%) d-|-Au events 

■\/^N]v = 200 GeV from the PHENIX experiment [2|. 
There is a reasonable agreement between the AMPT cal¬ 
culation and the V 2 data below px « 1.0 GeV/c, above 
which our calculation underestimates the experimental 
measurements. 

Leaving all AMPT parameters fixed, we show in Fig¬ 
ures and the V 2 and V 3 anisotropy moments as a 



Ab [rad] 
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FIG. 1. (a) Two-particle correlation for charged particles 

within 0.90 < px < 1.04 GeV/c for p+p and ^He-|-An events 
S'! \/^NN — 200 GeV. (b) Contributions to the correlation 
function arising from jet fragmentation are removed by sub¬ 
tracting away the per-trigger yield from p+p events. The 
resulting correlation function is shown in yellow. 


function of px for central p-|-Au and ^He-|-Au, respec¬ 
tively. The same general pattern of rising V 2 with px 
and smaller V 3 coefficients are observed for all systems. 
It is notable that there is an inflection point at px ~ 1-5 
GeV/c, after which both the V 2 and V 3 exhibit a slight 
decreasing trend. 

The RHIC experiments completed taking ^He-|-Au col¬ 
lision data at ^/s^ = 200 GeV in 2014. The PHENIX 
collaboration has presented V 2 and V 3 measurements in 
^He-I-Au at = 200 GeV |3]. These experimen¬ 

tal results are reproduced in Figure]^ and are found to 
be in reasonable agreement with the AMPT-extracted 
coefficients up to px ~ 1.0 GeV/c, beyond which the 
AMPT results fall significantly below the data, substan- 
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FIG. 2. AMPT calculation results for the azimuthal 
anisotropy moments V 2 and V 3 for central d+Au at 
— 200 GeV, compared with experimental results from the 
PHENIX experiment. 


tially more for V 3 than for V 2 if we consider the relative 
difference. 


IV. DISCUSSION 

Having obtained V 2 and V 3 for collision systems with 
different initial geometry and having compared them to 
experimental data, the question is then how to under¬ 
stand these results physically. In hydrodynamic models 
there is a straightforward physical picture of momentum 
anisotropy in terms of the velocity field of an expanding 
fluid. However, the situation is less clear in the case of 
transport models, such as AMPT. It has recently been 
proposed that the origin of the substantial V 2 calculated 
with AMPT lies predominantly in the anisotropic prob¬ 
ability of partons to escape from the partonic scatter¬ 
ing phase; that is, there is a preferential direction along 
which to freeze out M- In this section, we compare 
our elliptical and triangular anisotropy moments to un¬ 
derstand how they relate to intrinsic initial geometry in 
AMPT. 

In hydrodynamic modeling, the final momentum 
anisotropies are directly related to the initial state ec¬ 
centricities. In order to explore whether this relationship 
holds in AMPT, we take the ratios of AMPT calcu¬ 
lated V 2 and V 3 values between systems. Figure (top 
panel) shows the ratio of V 2 values as a function of pt in 
p-l-Au and ^He-|-Au relative to d-l-Au. Figure (bottom 
panel) shows the ratio of V 3 values as a function of px 



FIG. 3. AMPT calculation results for the azimuthal 
anisotropy moments V 2 and V 3 for central p-|-Au at 
= 200 GeV. 


in p-|-Au and ^He-|-Au relative to c?-|-Au. Also shown are 
the ratios of initial eccentricities, 62 and £3 from Table [ij 
as dashed lines. The AMPT V 2 and V 3 ratios between 
different systems are relatively px independent, with no¬ 
table deviations at lower px < 0.6 GeV/c, and following 
the same ordering as the initial geometry ratios. The 
values for the £ ratios are all lower than the AMPT 
ratios by approximately 15-35%. It is notable that the 
initial spatial eccentricities are calculated from the nu¬ 
cleons, and there will be minor variations in the geome¬ 
try upon string melting in AMPT. However, reasonable 
variations in the smearing parameter cr = 0.4 fm cannot 
resolve the full differences. 

The V 2 in p-|-Au is predicted to be substantially smaller 
than in d-l-Au and ^He-|-Au and the V 3 in ^He-|-Au is 
predicted to be substantially larger than in p-|-Au and 
d-l-Au, both from simple geometry. Although we have 
qualitatively confirmed the expected scaling of momen¬ 
tum anisotropy with initial geometry across a variety of 
small collision systems, the exact mechanism responsible 
for the anisotropies and the deviations from geometry 
scaling, for example at low px^ must lie within the inner 
workings of the AMPT model. 

Several authors have identified partonic scattering as 
the critical mechanism for the development of Vn in this 
context [HIIllIIH]. However, partonic scattering is only 
the first stage in the evolution of the partons that emerge 
from string melting, and it is relevant to examine the ef¬ 
fects that hadronization via coalescence and the subse¬ 
quent hadron cascade have on the final values. 

To that end, we repeat our calculations for ^He-|-Au 











5 



FIG. 4. AMPT calculation results for the azimuthal 
anisotropy moments V 2 and V 3 for central ^He+Au at 
— 200 GeV, compared with experimental results from the 
PHENIX experiment. 


and p+Au, enabling and disabling the hadronic cascade, 
and effectively turning off partonic scattering by setting 
the parton cross-section to very nearly zero. The results 
are shown in Figure]^ The first noteworthy feature in the 
top panels, corresponding to ^He-|-Au, is that keeping the 
partonic scattering phase with a = 1.5 mb, but turning 
off the hadron cascade has a substantial effect on In 
fact, the blue and violet curves show that a late-stage 
hadron cascade actually increases V 2 by about 20% for 
Pt > I GeV/c and by roughly 100% forpr < 0.5 GeV/c. 
This effect is even more pronounced for V 3 . Additionally, 
the green curves show that even when the partonic phase 
is turned off, a sizable V 2 and V 3 still develop by virtue 
of final-state hadronic interactions, with the effect being 
more pronounced for V 3 . Finally, as a consistency check, 
disabling both the parton and the hadron cascades results 
in the collapse of Vm as shown by the orange curves. 

The same analysis is carried out for p-|-Au in the bot¬ 
tom panels of Figure We observe that including a 
hadron cascade after the partonic scattering phase has 
a much smaller effect on the measured V 2 or V 3 . Fur¬ 
thermore, the Vn that develops from the hadron cascade 
alone when turning off partonic scattering is much less 
substantial than in the case of ^He-|-Au, as evidenced by 
the green curves, showing that in p-|-Au the bulk of the 
Vn originates in the parton cascade. 

We now examine the role of hadronization in the de¬ 
velopment of Vn- Figureshows V 2 and V 3 calculated for 
^He-|-Au and p+Au collisions with a partonic scattering 
phase using (z) partons at freeze out and (ii) hadrons 
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FIG. 5. (a) Ratio of elliptic and (b) triangular anisotropy 
moments as a function of pr in p-|-An and ®He-|-An compared 
with a baseline in d-|-Au. Dashed lines are the ratios of e 
values from the initial geometry. 


immediately after coalescence. For high pp, we observe 
that hadronization increases V 2 and V 3 in both collision 
systems. However, for pp < 0.5 GeV/c, the effect of 
hadronization is to reduce t;„, as evidenced by the cross¬ 
ing of the curves in the figure. This can be understood 
in terms of quark coalescence dynamics. Since hadrons 
are produced by aggregating partons in spatial proxim¬ 
ity and with collimated momenta, the coalescence yields 
hadrons with transverse momentum greater than that of 
their constituent quarks, hence increasing at higher 
pp. It is notable that this effect is greater in p-|-Au than 
in ^He-I-Au. 

The detailed mechanism and its relation to the num¬ 
ber of parton scatterings in the early AMPT stage still 
require further elucidation. That said, it is clear that the 
AMPT coalescence prescription and following hadronic 
cascade substantially modify and amplify this early stage 
effect. Therefore, the deviations from geometric scaling 
shown in Figure [^appear to arise from the relative dom- 
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FIG. 6. Effects on Vn of enabling and disabling the parton scattering phase and the hadronic cascade in AMPT for ®He+Au— 
panels (a) and (b), and p+Au—panels (c) and (d). 


inance of these stages as a function of px and collision 
system. 


V. SUMMARY 

Recent intriguing experimental observations at RHIC 
and the LHC have raised the question of whether small 
droplets of QGP can be formed in small collision sys¬ 
tems. From among several competing models, nearly in- 
viscid hydrodynamic calculations, both at RHIC and the 
LHC, give reasonable account of the measured anisotropy 
coefficients. However, parton scattering in transport 
models—AMPT with string melting, in particular—has 
also been shown to provide an adequate description of 
the long-range azimuthal correlations and momentum 
anisotropy coefficients measured in p-\-p and p-|-Pb at 


LHC energies. In this paper, we extend these calcula¬ 
tions to RHIC energies, focusing on the insight that can 
be gained by varying the initial geometry of the projectile 
nucleus. 

We find that AMPT is capable of reasonably repro¬ 
ducing the measured elliptic and triangular flow coeffi¬ 
cients for central d-|-Au and ^He-|-Au collisions at 
= 200 GeV for pp < 1 GeV/c. However, AMPT un¬ 
derestimates the measured values for higher px- With 
this observation, we ascertain the validity of the model 
for rendering initial geometric anisotropy into final-state 
particle momentum correlations for small systems at both 
the RHIC and LHC energy scales. We also make predic¬ 
tions for elliptic and triangular anisotropy coefficients in 
p-l-Au collisions at = 200 GeV and qualitatively 

relate these results to calculated initial-state geometric 
anisotropy. 
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FIG. 7. Elliptic and triangular anisotropy moments in (a) ^He+Au and (b) p+Au collisions following partonic scattering using 
partons immediately prior to, and hadrons immediately after hadronization via coalescence. 


However, we also find that partonic scattering is not 
the only source of the substantial elliptic and trian¬ 
gular momentum anisotropies in the AMPT model. 
Hadronization and the subsequent hadron cascade exert 
important modifications on the from partonic scat¬ 
tering, with strong dependences both on pT and the in¬ 
trinsic initial geometry of the system. Direct compar¬ 
isons with experimental data in these new systems, with 
both AMPT and various hydrodynamic models, is antic¬ 
ipated to shed light on the physical dynamics involved in 
these collision systems. To finalize, we highlight the need 
to identify additional observables that provide a more 


stringent discrimination between initial geometry and its 
translation to final-state correlations. 
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